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ABSTRACT 

The  probabilities  of  single  and  two  successive  gas  surface 
collisions  satisfying  minimum  interaction  energy  cutoffs  are  considered 
in  the  hard  cube  model  sense.  Comparison  is  made  with  atomic  and  mole- 
cular oxygen  reaction  probabilities  with  graphite  corresponding  to 
specific  mass  ratio  examples  of  4/3  and  8/3.  Comparable  peak  probabili- 
ties result  for  multiple  interactions  but  of  much  reduced  sharpness  with 
temperature  for  only  a two  collision  constraint.  Realistic  activation 
energies  are  implied  from  a match  of  the  energy  cutoff  probability  model 
to  the  portion  of  the  experimental  data  increasing  with  surface  tempera- 
ture. The  rapidly  decreasing  reaction  probabilities  found  experimentally 
at  higher  temperatures  imply  an  appreciable  increase  in  the  required 
number  of  effective  collisions. 
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SYMBOLS 


A, 

A,  B,  C,  D 
G 

H 


‘j  <%> 

Oj  (^o) 

k 

k.. 


4 


"j 


m* 

N 

dP 


Toe  X 


Pl(F) 

Pjj(e^.E2) 


regions  of  (Uq,  Vq)  plane  used  for  Pjj  evaluations 
integration  limit  parameters  for  I.,  Equations  (40),  (41) 

J 

distribution  function  for  surface  particle  speed  |Vq|, 
Equation  (1) 

distribution  function  for  gas  particle  normal  velocity 
component  u^.  Equation  (2) 

integral  form.  Equations  (40),  (42) 

integral  form.  Equation  (43) 

Equations  (45),  (48) 

Boltzman  constant 

slope  of  border  to  B^  region 

Equation  (58) 

surface  "box"  length;  also  integration  limit  for  Jj 

Equations  (50),  (51) 

gas  particle  mass 

surface  particle  mass 

reduced  mass,  m^  m^  /(m^  + m^) 

Equation  (49) 

surface  particle  location  probability.  Equation  (20) 

directional  probability  that  first  collision  is  inward 
or  outward  for  given  (u^,  v^)  pair.  Equation  (25) 

probability  that  second  collision  follows  from  (u  , v ) 
first  collision 

probability  o£  first  collision  with  interaction  energy 
in  excess  of  Ci 

probability^ of  second  collision  with  interaction  energy  in 
excess  of  ^2  for  first,  second  collisions  respectively 
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Pi,  Qi.  Ri 


T 

T* 

U 

u„ 


Ui 


‘^0’  "o 


0 coll 


01 


X* 

(xVL) 


m 


y 


coefficients  in  integrand  of  region  contribution  to 
Pjj,  Equation  (65) 


incident  gas  temperature 
surface  temperature 


Ts/Tg  cos  0,  Equation  (41) 


(mT) 


integration  limit  for  J. 


gas  particle  normal  velocity  component  prior  to  first 
collision,  positive  inward 


gas  particle  normal  velocity  component  after  first 
collision,  positive  inward 


. . 1/2  1/2 

normalized  particle  velocities,  Xg  u^,  x^  v^. 


surface  ,'article  velocity  prior  to  both 

first  collision  and/or  rebound,  positive  outward 


surface  particle  velocity  immediately  prior  to  first 
collision,  positive  outward 


Vq  on  cut-off  energy,  e^- , locus.  Equations  (8),  (17) 


surface  particle  velocity  after  first  encounter, 
positive  inward 


distance  from  outer  edge  of  surface,  positive  inward 
first  collision  location 


minimum  location  (normalized)  for  first  collision  which 
allows  second  collision  possibility.  Equation  (27) 


m„/2kT„cos  0 

g g 


m3/2kl3 


Equation  (52) 

Cut-off  energy  for  first  collision 
nondimensional  cut-off  energy.  Equation  (11) 
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Et,  e,»  first,  second  collision  relative  integration  energy, 

^ Equations  (7),  (12) 

(cp).  reaction  probability  for  oxygen  on  graphite,  j = 0,  O2, 

^ ^ Equation  (72) 

e*  [(1  + 

0 incidence  angle  measured  from  surface  normal 

y ratio  of  gas  to  surface  particle  masses,  m /m 

9 ^ 
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CHAPTER  1 
INTRODUCTION 


Measured  reaction  probabilities  for  both  molecular  and  atomic 

12  3 4 

oxygen  reactions  with  carbon  show  anomalous  (peaked)  Arrhenius  plots.  ’ 

Comparison  of  results  for  the  molecular  and  atomic  oxidation  cases  indi- 

1 4 

cate  the  atom  reactivity  to  be  higher  by  an  order  of  magnitude.  ’ It 

has  been  suggested  that  the  increase  may  be  related  to  a higher  sticking 

probability,^  and  that  dissociative  adsorption  may  be  the  limiting  step 

2 4 

in  the  reaction  mechanism.  ’ 

There  is  some  evidence  that  primarily  normal  energy  barriers 
are  important  in  the  dissociative  adsorption  and/or  recombinative  desorp- 
tion processes. In  simplest  form,  interactions  with  greater  (normal) 
energies  than  some  cut-off  energy  will  then  lead  to  adsorption  in  dissoci- 
ated form.  For  a gas  at  rest  with  respect  to  an  adjacent  solid  surface 
the  intensity  of  all  particles  impinging  on  the  surface  varies  as  the 
cosine  of  the  angle  from  the  surface  normal.  However,  the  intensity  is 
non-cosine  for  those  particles  with  normal  energies  greater  than  some 
non-zero  cut-off,  and  the  distribution  is  then  peaked  in  the  normal 
direction. 
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These  observations  add  to  the  importance  of  the  normal  inter- 
action, which  if  dominant  offers  considerable  simplification  in  modeling 
the  events  of  gas  and  solid  particle  encounters.  The  hard-cube  model 
has  had  some  success  with  precisely  that  assumption  assuming  surface 
particles  confined  to  a rigid  box  and  interacting  with  the  normal  compo- 
nent of  Maxwellian  gas  particles.  However,  a single  collision  restriction 
was  imposed  by  limiting  the  mass  ratio  of  gas  and  surface  particles  to  be 
p < 1/3.  Primary  interest  at  the  time  was  in  inert  gas  scattering 
distributions,  for  which  there  was  comparison  data  from  experiment  and 
with  which  there  proved  to  be  quite  reasonable  agreement. 

For  reaction  encounters  the  natural  extension  of  a hard-cube 
approach  is  an  allowance  for  an  energy  barrier  and  removal  of  the  p < 1/3 
constraint.  The  former  corresponds  to  imposing  a reaction  probability, 
while  the  latter  introduces  multiple  collisions  which  may  weight  the 
probability  according  to  changes  in  the  residence  time  on  the  surface. 
Single  and  double  collision  interactions  are  considered  here  while 
retaining  the  simplifying  assumption  of  negligible  lateral  interaction. 
Thus,  near  normal  encounters  are  implied  and  with  only  one  atom  of  the 
solid  surface. 
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CHAPTER  2 

COLLISION  PROBABILITIES  MODEL 

2.1  Distribution  Function  Basis  for  Probability 

Goodman^^  points  out  that  for  p = nig/m^  < 1/3  certain  types 
of  two  successive  collision  sequences  are  ruled  out,  and  that  for 
p < 1 there  will  be  two  collisions  at  most.  However,  an  oxygen-graphite 
encounter  corresponds  to  p > 1 , (i.e.,  either  16/12  or  32/12)  and  multiple 
collisions  must  be  expected  here. 

We  assume  a constant  speed  surface  particle  moving  normal  to 
the  surface  and  instantly  reversing  direction  at  the  ends  of  its  one 
dimensional  box  enclosure. Immediately  prior  to  a first  collision 
its  velocity  is  v^  positive  in  the  outward  direction,  and  of 

magnitude  0 < jv^  ^ The  probability  of  a first  collision 

involving  |v^  ("IvqI)  is  the  probability  that  the  speed  is  |Vq|  at 

the  time  the  gas  particle  enters  the  box.  From  the  one  dimensional 
Maxwellian  distribution  this  is 

G(Ivj)^k  = 0) 

where  ^ • 

An  approaching  gas  particle  in  a Maxwellian  beam  oriented  at 
an  angle  0 from  the  surface  normal  has  a normal  velocity  component 
positive  inward,  with  the  distribution  (o  s u^  i «) : 


(2) 
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r 


Classical  encounters  will  be  assumed  between  the  surface  and 
gas  particles  along  the  normal  direction.  Tangential  interactions  are 
neglected  and  clearly  this  is  a progressively  poorer  approximation  with 
an  increasing  number  of  collisions  prior  to  departure  of  the  gas  particle 
from  the  box.  Indeed,  for  a sufficiently  large  number  of  collisions  it 
should  be  expected  that  a substantial  part  of  the  tangential  momentum  is 
absorbed  for  any  one  gas  particle.  Should  a second  collision  occur  it  is 
assumed  to  be  with  the  same  surface  particle.  The  tangential  and  second 
collision  assumptions  together  effectively  limit  the  range  of  0 to  "near" 
zero.  The  explicit  ap  earance  of  0 is  then  solely  in  x^,  and  in  the 

P 

combination  T^cos  0,  and  is  effectively  then  a beam  cooling  parameter. 

Consider  the  probability  that  a first  collision  occurs  for  a 
specific  interaction  energy.  For  such  collisions  to  occur  prior  to 
rebound  of  the  surfact  particle 

either,  'Vo  ^ ^ IL  > , for  an  outward  type 

cea  ^ 

or,  ^ < O , for  an  inward  type  (3) 

whereas  if  rebound  precedes  the  collision, 

'Vq^  - ~Vo  . for  an  outward  type  (4) 

This  may  ‘ - interpreted  as  an  increased  directional  probability, 
PplUo,  Vg),  for  the  v^  > 0,  i.e.,  outward  type,  collisions  and  a corre- 
sponding decrease  in  the  probability  for  v^  < 0,  i.e.,  inward  type,  colli- 
sions. With  that  understanding,  the  probability  of  a first  collision 

0 coll  ” '^0 


where 


with  V 
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1 


dPj  = '==4'^..  (5) 

at  a point  in  the  (u^,  v^)  half  plane.  Only  some  fraction  of  these  will 
result  in  a second  collision.  If  P2^^o’  '^o^  probability  of  a 

second  collision  resulting  from  an  initial  (u^,  v^)  pair  for  the  first 
collision,  then 

^Pji:  ~ 

is  the  overall  second  collision  probability  for  a point  in  the  (u^,  v^) 
half  plane.  Integration  of  Equations  (5)  and  (6)  over  portions  of  the 
half  plane  (0  s u^  < <»;  ^ v^  ^ <»)  furnishes  the  proportional  number 

of  first  and  second  collisions  relative  to  the  number  of  gas  particles 
arriving  at  the  surface. 

In  view  of  the  Pq  weighting  in  Equation  (5)  only  Equation  (3) 

(i.e.  Vq  ^qii  = Vq)  need  be  considered  and  the  relative  interaction 

1 2 

energy  for  a first  collision  is 

e^.  =.  (7) 

in  terms  of  a reduced  mass  are*  = ^ '"g^’  '^o  ^ ”*^o 

all  lesser  speeds  correspond  to  Equation  (4).  For  relative  energies 
equal  or  greater  than  a cut-off  energy  c-j , say,  Equation  (7)  imposes 
the  constraint  (Fig.  1): 

The  relative  number  of  such  collisions  is,  then,  from  Equation  (5), 


(8) 
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ffO  K) 


O u. 


(9) 


It  is  worth  noting  the  limiting  case  of  T -*■  0 for  which  |v  | -*-  0 and 
Cir 


^-^2e,/A<*  ; i-e.. 


<30 


^2  f^). 


A X -1 

=A’'r^V'v/i? 


Z_ 


(10) 


in  which 


e-  = — 


(11) 


is  the  natural  nondimensional  interaction  cut-off  energy. 

A second  collision  energy  cut-off  follows  from  the  relative 
interaction  energy 


(12) 


where  u^ , are  the  particle  speeds  immediately  after  the  first  encounter. 
The  sign  choice  in  Eguation  (12)  follows  from  v-|  > 0 for  p > 1 since  from 
momentum  and  energy  considerations 


U/ 


(13) 
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if  a collision  occurs  (v^  ^ ^1’  '^1  positive  inward. 

As  a consequence,  the  second  collision  occurs,  if  at  all,  after  the  sur- 
face particle  rebounds  from  the  bottom  of  the  box,  and  the  relative  speed 
is  then  (u.j  + v-j)  for  u^j  ^ 0.  From  Equation  (13)  a second  collision 
always  (i.e.,  ~ takes  place  for 


V , J 


(14) 


and  takes  place  for 


V^>/^Uo 


■Y  <^c>  j ~ j (15) 

only  if  the  gas  particle  is  overtaken  prior  to  its  departure  from  the 
surface  (i .e.  ?2  < 1.0). 

A cut-off  constraint  for  int“*"t»ction  energies  ^ C2  is  then 


7^ 


and  corresponds  to  (Fig.  1): 


? - 


(16) 


y 


3-/ 


(W) 


for  u ^ 3 on  introducing  Equation  (13).  The  T -*•  0 limit  here  implies 


> 


/y 


(18) 


for  both  \i  i 2.  Since  P2  = 1 for  |Vq|  -*■  0 in  view  of  Equation  (14), 
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■) 

=y  /y'^4',  = [/i- e 


U-c 


where  \i-^  = ^ , ^2  ~ 3^] y » Equation  (19)  requires  a lower  limit  of 

integration  choice  which  is  the  larger  of  the  u^  minimums  [either  Equation 
(18)  or  that  above  Equation  (10)]. 


2.2  Directional  Probability,  Pr 


At  the  time  that  a gas  particle  enters  the  surface  box  the 
solid  particle  has  neither  preferred  direction  nor  location.  Thus  the 
probability  that  the  initial  motion  is  either  outward  (v^  > 0)  or 
inward  (v^  < 0)  is  1/2,  and  the  initial  location  within  a box  of  length 


_ /^/ 


However,  our  interest  is  in  the  directional  probability  at 
the  time  of  first  impact,  Pp,  such  that  the  collision  itself  is  of  an 
outward  (v^  > 0)  or  inward  (v^  < 0)  type  for  a given  (u^^,  v^) 

pair.  Clearly  all  > 0 will  result  in  outward  collisions.  In  addi- 

tion some  Vq  < 0 [see  Equation  (4)]  also  lead  to  outward  collisions  if 
rebound  occurs  prior  to  being  overtaken  by  the  gas  particle.  It  follows 
that  inward  collisions  will  occur  for  only  the  remaining  < 0 particles, 
If  -Uq  < Vq  < 0 the  times  required  for  the  surface  and  gas 
particles  to  reach  the  bottomeof  the  box  are  (L  - x)/|Vq|  and  L/u^ 
respectively.  The  initial  locations 
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<r  / _ 

4 > ‘‘o 

(21) 


therefore  imply  overtaking  (inward,  < 0)  and  rebound  (outward, 

Vo  (,qii  > 0)  collisions  respectively  when  v^  < 0.  Combining  the  initial 
direction  probabilities  with  Equation  (20)  for  the  box  regions  of 


Equation  (21): 

r 0 


('-a^<Vo<o)  (22) 

( O <Vo  } 


( o <iro  ) 


for  any  v^  the  probability  of  a first  collision  is  ^Din  ^ ^Dout  so  that 


and  indicates  Pj  to  be  unity  for  the  entire  half  plane  as  would  be  expected. 

With  energy  constraints  such  as  Equation  (8)  it  is  more  conven- 
ient to  collect  the  inward  and  outward  contributions  together  to  obtain 
the  equivalent  directional  probability  as  in  Equation  (5).  I.e., 
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to  indicate  outward  type  collisions  (v^  > 0)  are  more  probable  for 

a given  |Vq|,  and  that  the  overall  probability  for  any  |Vq|  is  in  fact 
unity  for  a first  collision  (Figure  2).  E.g,,  for  the  entire  half  plane 

Equation  (24)  may  then  be  written  as 


2.3  Second  Collision  Probability,  Pq 

A second  collision  is  certain  if  v^  < (u  - 1 )/2Uq  in  view 
of  Equation  (14).  Thus 

i /t-W  (26) 

The  remaining  portion  of  the  (u^,  v^)  half  plane  involves  only  outward 
Vq  > 0 first  collisions  with  u-j  < 0 as  in  Equation  (15).  Some  of  these 
collisions  take  place  at  a location,  x*  say,  which  permits  the  surface 
particle  to  then  rebound  from  the  bottom  of  its  box  and  overtake  the  gas 
particle  for  a second  collision.  P2  follows  from  the  probability  that 


L 


a first  collision  will  occur  within  the  allowable  x*  range  consistent 
with  the  arbitrary  surface  particle  location.  Equation  (20),  and 
direction  initially. 

Allowable  x*  locations  are  those  for  which  the  gas  particle 
exit  time,  x*/(-u-|),  exceeds  the  equivalent  time  for  the  surface  parti- 
cle, (2L  - x*)/vy  Using  Equation  (13),  overtaking  second  collisions 
correspond  to  first  collisions  in  the  range 


(27) 


Larger  (v^  (•oll'^^^  imply  a smaller  range  of  allowable  locations.  For 
(Vo  coll'^'^o^  = (y  - D/2  the  collision  may  occur  anywhere  in  the  box, 
consistent  with  Equation  (26).  For  \i  > 3 there  is  an  allowable  range 
for  all  (Vjj  coll''%'  < ^or  p < 3 a sufficiently  large  v^  coll^^o 
[>  (3y  - l)/(3  - u)]  cannot  result  in  an  overtaking  second  collision, 
consistent  with  the  energy  constraint  in  Equation  (17)  when  C2  = 0. 

Since  outward  (v^  > 0)  first  collisions  arise  from  both 

Vq  > 0 the  collision  site  differs  for  each  initial  direction.  Equating 
the  times  for  the  two  particles,  the  first  collision  occurs  at 


'^0>0  ) 


(28) 
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The  probability  that  the  collision  occurs  within  d(x*/L)  is  then 


JP  -=1^ 


h^\  ^ 

A.  / 4 


(Voio) 


(30) 


and  the  collisions  are  in  the  ranges,  from  Equation  (29), 


The  limitation  on  initial  locations  x/L  for  the  inward  directed  surface 
particles  with  velocities  between  0 > v^/u^  > -1  corresponds  to  those 
particles  that  rebound  and  thereafter  undergo  an  outward  (v^  > 0) 

interaction  (Equation  21),  Note  that  there  is  an  implicit  influence  of 
the  mass  ratio,  p , in  Equation  (31)  by  virtue  of  the  present  interest 

(''o  coll^^)  " 

The  second  collision  probability,  P^,  follows  from  Equation  (30) 
applied  to  those  portions  of  the  first  collision  ranges,  Equation  (31), 
which  fall  within  the  allowable  range  for  a second  collision.  Equation  (27). 

For  the  direct  outward  (v^  > 0)  interactions:  [(p  - 1 )/2  < v^/Up] 
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f-i  - o 


For  the  rebound  outward  (v^  < 0)  interactions;  for  (y  - l)/2  < 1 


If 


(32a) 


(32b) 

< lVol/S  = 


(33a) 


(33b) 


(33c) 
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and  for  (y  - 1 )/2  < |Vq|/Uq  < 1 : 


(34a) 


(34b) 


2.4  Overall  Second  Collision  Probability,  Pjj 

The  mass  ratio  y provides  the  basis  for  specific  integration 
intervals  when  evaluating  Pjj  from  Equation  (6)  with  (25).  Note  that 


■ >/<  > 

> 

(35) 

¥ 

$ / 

s 

(36) 

2 

^ / 

3 
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Then  for  1 < y < 5/3,  5/3  < y < 3,  and  3 < y respectively: 


prior  to  any  energy  cut-off  considerations.  The  essential  point  is  the 
association  of  a proper  accord  with  the  latter's 
direct  or  rebound  first  collision  origin. 

For  the  specific  interaction  of  atomic  oxygen  and  carbon,  i.e.. 


y = 16/12,  the  contributions  from  Equations  (25),  (26)  and  (32)  - (34) 
into  (37)  gives  (for  e.|  = C2  = 0): 


Similarly  for  molecular  oxygen,  y = 32/12: 


where  for  each  interval  for  which  > Ug(p  - 1 )/2  the  explicit 
'’o^direct  ^ (^2  '"oVebound^  indicated. 

The  energy  constraints  v„  > v„  and  v„  < affect  only  the 

O 0 1 0 O2 

integration  limits  as  in  Figure  1.  Explicit  integral  results  are  deve 
loped  in  the  next  section  in  appropriate  nondimensional  form. 
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28 


(44) 


(45) 


(46) 

(47) 


(48) 


(49) 
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and  the 


(50) 


are 


M_  = 


' ■ - vS.^'h? 


Mo  = 


tt^h 


Zh^f 


(51) 


y 


/ cf 
O-fc^h 


Here 


cf 


(t+c’'js^  / CP 


C/^c‘)'^ 


(52) 


The  special  limit  of  (T5|  -*■  <»  then  corresponds  to 


M„  = M,  = M,  = K,  = 0 
0 1 2 I 


- I 


•^0*  7 


(53) 


K,  = 


V ^ 


— L 


+ T-2-e 


-u. 


1 
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and  application  to  Equation  (24)  yields 

= ^ [o  - -^  / -&/J  ~ I (54) 

as  a simple  example. 

3.2  First  Collision  Probability,  c-j  > 0 

For  a finite  energy  cut-off  as  in  Equation  (8)  the  first 


collision  probability  is  (compare  Equation  25a) 


Equation  (55)  in  nondimensional  form  analogous  to  Equation  (42)  then 

rpHiiroc  tr> 
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where 


(58) 


1/^.  = (59) 

Equivalently,  for  the  three  regions  a,  b,  and  c as  shown  in  Figure  3: 

-u> ] 

(Ah 

(Ah  ‘ [w 


and  the  Pj  (e*,  T*)  dependence  may  be  interpreted  as  Pj  (p,  ■£ -j , T)  on  the 
basis  of  Equation  (56).  Thus  the  corresponding  energy  and  temperature 
ratio  pairs  for  specific  collision  probability  are 

~ 

for  distinct  particles  and  P|^ 


3.3  Overall  Second  Collision  Probability,  e.|  > 0 C2  > 0 

p 

II  follows  from  Equation  (37)  modified  for  finite  cut-off 
energy  levels.  For  example,  for  p < 3,  i.e.  for  dv„  /du„  > 0 for  the  e, 

02  O 

cut-off  locus  defined  by  Equation  (17),  two  types  of  integration  regions 
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provide  contributions  to  Pj-j.  These  will  be  referred  to  as  Type  A and 
regions,  and  in  each  case  the  result  will  be  seen  to  depend  upon  the 
relative  cut-off  energy  ratio,  ca/ci,  that  implies  limitations  on  the 
relevant  portions  of  the  regions  (Figure  4).  Region  boundaries  consist 
of  the  v„  and  v„  cut-off  locii,  and  the  natural  boundaries  implied  by 

Oi  O2 

as  in  Equation  (37),  i.e.,  v^  = k^u^  with  = (y  - l)/2,  (3y  - l)/4, 
1 and  y. 

The  contribution  from  Region  A,  for  which  k^  = (y  - l)/2 
depends  upon  the  relative  positions 


or,  in  general  for  any  k.. 
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The  contribution  from  a Region  requires  consideration  of 
(^2/ei)  relative  to  both  {c2/ei  and  (e2/ei)^  + 1 -phus 


Si  7fi^ 


C^g/  ) 


"‘’I 


J 


J'<^^(65a) 


and  depending  upon  (^o^c^^^o^bi 

Aa\  N^Ju 


^^IcSCsiOU  6 c 

,^l4i 

(W  ^ ^ 


ii 


^ a,t1  ^ 


S<^%dUo 


xn; 

lUi 


(65b) 
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where 


5 = 
B 


z 


lU  t 


4 


f/^}  [ 

Also,  U and  7 , are  given  by  Equations  (62)  and  (59),  and  the  remaining 


oa 

limits  are: 


Ifo. 


Since 


f^-/)  - 

(/■hu) 


(66) 


it  follows  that  for  (02/^1)  > 1 only  Equations  (64b)  and  (65d)  need  be 
considered.  The  minimum  energy  criterion  entering  into  Equation  (65)  is 
that  corresponding  to  the  = yu^  upper  bound  (see  Equations  (38)  and 
(39) , e.g. ) for  which 

ik]  ^ 


(67) 


(68) 


from  Equation  (63).  For  lesser  energy  ratios  only  Equations  (64a)  and 
(65a)  need  be  considered  and  the  resulting  Pjj  f Pjj(c2). 

For  the  special  case  of  (cz/ci)  > 1 Equations  (64b)  and  (65d) 


reduce  to: 


[ I --J/fej  ^ -Z/Kjj 


-r*  rV 


f.i,  - 


-r",TV 


"b  1 


*f  ^c  [_^>(^c)~-i'i(^'] "f" ^ [^(^^)~X^)Jj  . . 

I ^ 

f [xw -r/ife)]+ ^ frhl -J  1 


From  Equations  (38),  (39),  and  (63),  for  the  specific  mass  ratios  of 
interest  here: 


U = 4/3 


U •=  8/3 


'8 

i 

1 

I 2 6 

4 

" 

44  , 

0 

0.0816 

: 1 

1 

0 

-3 

8 !o.207 

1 

11 

ITT 

4 ' 

0.0204 

1 
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In  general,  Pjj  = Pjj  (e*>  T*,  y,  e^/e-i),  and 

” t i 

The  results  for  Pj,  Equation  (60),  and  Pjj,  Equations  (64), 
(65),  were  progranned  for  numerical  evaluation  over  a range  of  e^-  and  T 
for  y = 4/3  and  8/3.  The  computer  listings  are  provided  in  Appendix  1. 
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CHAPTER  4 
DISCUSSION 

4.1  Reaction  Rates  and  Energy  Cut-offs. 

4 

A reaction  probability  may  be  defined  as 

carbon  atom  flux  from  surface  /J2) 

^R  ” incident  oxygen  flux  to  surface  ' ' 

and  implies  maximum  (Ep)Q,  values  equal  to  1 and  2 respectively 

for  atomic  and  molecular  beams  which  result  in  CO  as  the  reaction  pro- 
duct. The  results  of  Liu's^  measurements  are  given  in  Table  1 and  (e[^)g 
and  (e[^)Q^/2  may  be  compared  with  the  collision  interaction  energy  prob- 
abilities outlined  in  Chapter  3 as  a measure  of  the  energy  constraint 
importance  for  effective  interactions. 

Liu  used  an  oxygen  beam  at  9 = 45°  incidence  on  a graphite 
target  at  temperatures  in  the  range  900  ^ T^,  °K  < 1800.  From  beam 
temperatures  of  T^  = 300,  1475,  1875,  and  2200°  K,  he  inferred  the  sepa- 
rate atomic  and  molecular  reaction  probabilities,  {^^)q  and  (e(^)g^,  the 
dissociated  case  being  based  on  the  highest  beam  temperature  (2200°  K) 
for  which  the  Knudsen  source  provided  a 15%  equilibrium  dissociation  level. 
Table  1 indicates  (ep)g^  to  have  a primary  dependence  on  surface  tempera- 
ture, T^,  for  the  appreciable  gas  temperature,  T^,  range  considered. 

The  probability  of  an  encounter  with  minimum  interaction  energies 
CpCj,  is  of  the  form  P(e^,  T,  y,  with  the  normalized  energy  re- 

ferred to  Tg  as  in  Equation  (11). 

For  0 = 45°  this  implies 
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e.  = (0.993  X 10'^)(Tg,  °K)e.  * (lO'^Tg)e^  kcal/mole  (73) 

and  therefore  provides  some  indication  of  cutoff  energy  level,  , 
variations  with  Tg  if  constant  are  indeed  governing  for  any  collision. 

4.2  Collision  Probability  Results. 

The  fundamental  first  collision  probability,  Pj(E].  T,  y)  is 
shown  in  Fig.  5 for  y = 4/3.  Tabulations  of  the  calculated  results  for 
all  figures  are  included  in  Tables  2 through  5.  The  probability  of 
achieving  a given  energy  cutoff  decreases  with  increasing  e.j , and  de- 
creasing T as  follows  from  the  implicit  (decreasing  T^,  fixed  T^)  and 
(decreasing  T^,  fixed  Tg)  bases  respectively.  Of  interest  first  is 
the  resulting  range  1.0  > Pj  > 0.1  corresponding  to  ^ 5,  or  e-j  'v 
0(10  kcal/mole)  from  Equation  (73).  Secondly,  the  probability  variation 
near  T-v  0(1)  proves  to  be  quite  appreciable  but  provides  no  basis  for 
a peak  reaction  level  as  is  evident  from  the  measurements. 

However,  increasing  T implies  relatively  larger  v^  magnitudes 
for  all  interactions,  and  Equation  (15)  then  implies  that  a lesser  number 
of  subsequent  collisions  may  then  be  attainable.  I.e.,  the  larger  nega- 
tive velocity  imparted  to  gas  particles  results  in  a larger  number  of 
escapes  prior  to  being  overtaken  by  a rebounding  surface  particle.  Fig. 

6 indicates  the  effect  on  Pjj  of  a minimum  combined  constraint  for  two 
collisions,  (e  ^ E2)  = (0,  0),  and  suggests  a resulting  basis  for  the 
appearance  of  a peak  probability  if  a multiple  interaction  is  in  fact  a 

_3 

necessity  to  affect  a reaction.  Similarly,  the  (1,0)  and  (10  , 1)  constant 

locii  indicate  the  scale  of  the  interaction  to  be  as  small  as  0(1) 
kcal/mole  in  such  cases  if  P > 0.1. 
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Figures  7 and  8 show  the  peak  development  explicitly.  For  small 
T,  F.|  dominates,  whereas  for  large  T,  £2  dominates.  For  the  two  collision 
model  peaking  occurs  for  T-v  0(1)  for  energy  levels  up  to  0(10) 
kcal/mole.  Clearly  an  increasing  number  of  collisions  implies  sharper 
peaks  and  their  occurrence  at  lesser  T. 

Figures  9 and  10  indicate  similar  Pj  and  Pjj  probabilities  for 
y = 8/3.  The  lesser  (£|^)g^  values  result  in  larger  choices  but  quali- 
tatively similar  behavior.  Larger  mass  ratio,  however,  implies  a shift  of 
peak  probabilities  to  relatively  larger  T levels,  as  must  be  expected 
from  the  increasing  demands  on  Vq  in  the  inequality  Equation  (15). 

4.3  Energy  Levels. 

Some  activation  energies  reported  for  the  atomic  oxidation  of 
14 


graphite  are  as  follows 

to  900 
473  - 573 
488  - 573 
287  - 473 
293  - 573 


kcal/mole 

0 

0-10 

6.5 

10 

1-13 


.15 


For  molecular  oxidation  the  comparable  activation  energy  is 
in  the  range  55-65  kcal/mole. 

Figure  5 shows  a first  collision,  constant  probability 
variation  with  T that  is  consistent  with  the  increasing  experimental 
reactivity  up  to  T %1300°  K. , and  Figures  6-8  indicate  marked  reductions 
in  probability  for  any  accompanying  second  collision.  The  fit  with  expe- 
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riment  of  the  calculated  e-j  % 2.5  locus  corresponds  to  5.5  kcal/mole  and 

14 

is  in  reasonable  agreement  with  the  above  reported  energies  for  p = 4/3. 

The  several  data  groups  obtained  for  distinct  beams  in  the 
p = 8/3  molecular  case  imply  first  collision  energy  level  requirements  of 
magnitude  (Fig.  10): 


e-j  , kcal/mole 

9. 

13. 

13. 

The  approximately  constant  e-j  values  above  confirm  the  expected  shift  of 
experimental  data  with  T.  This  appears  to  add  some  importance  to  the 
first  collision  energy  criteria. 

For  both  atomic  and  molecular  oxidation  it  is  evident  that 
multiple  collisions  are  required  for  consistency  with  a decreased  reacti- 
vity at  higher  surface  temperatures.  This  is  not  totally  inconsistent 

2 1 g 

as  well,  with  suggestions  of  surface  annealing  and  activation  ’ as 
competing  processes  that  tend  to  dominate  the  reaction  control  at  high 
and  low  surface  temperatures  respectively.  A less  reactive  surface  at 
higher  temperatures  necessarily  should  require  either  a greater  first 
collision  energy  cut-off  or  represent  a longer  dwell  time  on  the  surface 
during  which  further  collisions  can  take  place.  These  results  indicate 
a marked  decrease  in  probability  must  be  expected  if  a second  collision 
is  required.  For  example,  for  p = 4/3  and  near  the  peak  reaction 
probability  (T^  ^ 1300°  K)  temperature,  Pj(2.5)  !\s  2.2  and  Pjj(2.5,  0) 
0.065. 


T 

300 

30 

1475 

9 

1875 

7 

(. 
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4.4  Conclusion. 

The  probability  of  encounters  between  gas  and  surface  particles 
of  sufficient  energy  provides  some  of  the  qualitative  features  found  for 
oxygen-graphite  measured  reaction  probabilities.  Specifically,  increased 
reactivity  with  increasing  surface  temperature  may  be  attributed  to 
constant  interaction  energy  constraints,  and  sharply  decreasing  reacti- 
vity at  higher  temperatures  may  be  associated  with  the  very  much  reduced 
probabilities  that  follow  if  any  additional  collisions  are  necessary. 


f 
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APPENDIX  1 
COMPUTER  LISTINGS 

Listings  for  the  individual  programs  for  one  and  two  collision  proba- 
bilities are  given  below.  In  each  case  only  a single  input  data  card  is 
shown  (preceding  the  last  deck  card,  "/*")  as  an  example  of  the  multiple 
input  cards  called  for  by  the  J = 1 , 4 or  J = 1 , 26  counters. 

The  symbol  equivalence  is: 


EB 

= 

^1 

TS 

= 

I* 

EBR 

= 

UA 

= 

EBRI 

= 

( £3/^1 ) 1 

UC 

= 

("o’c 

EBRIl 

= 

+ 1 

UE 

= 

<%>e 

ES 

= 

e* 

V2I 

= 

PI 

= 

V2R 

= 

slope  of  v„  (IT 

O2  c 

P2A 

= 

^^II^A 

XMU 

= 

p 

P2BB 

= 

^11 

XKI(j) 

= 

•'i 

P2B(I) 

= 

(’’lI^Bi 

XMj 

= 

Mj  (j  = 0,1,2) 

PI 

= 

Pi 

+ 

X 

= 

Kj  (k  = 0,1,2) 

QI 

= 

Qi 

Y 

= 

y 

RI 

= 

Pi 

YN 

= 

N 

TB 

= 

T 

The  two  collision  listing  specifically  provides  for  either  of  two  XMU 
(=  p)  by  including  the  appropriate  constants  [see  table  after  Equation 
(70)]  in  cards  labeled  20  to  28. 


* i-TiKCKDI.G  iUTi  ^LA!lK-ncT 
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c 

c 

c 


c 


c 


c 


ONE  COLLISION,  BASIC  TEST.  FINITE  0,0 


dimension  XK{3),  XP(8) 

common  XK,  TS 

INPUT 

00  53  J=  1,A 

4 READ  5,  XMl),EB,TB 

5 EORMAT  (FIn.7,?F7.3) 

ES  = SORT  ((1.  ♦ XMU)  • 
ES2  = ES/2. 

TS  = l./SQPT  (XMU  • TB) 
ET  = ES  • TS 


7 CALL  ETEL 
9 CALL  EYEL 
call  EYEL 
CALL  EYEL 
CALL  EYEL 
CALL  EYEL 
9 call  EYEL 
10  CALL  EYEL 


(-TS,ET.O.,B.,ESa,l,xP(n I 
(-TS,ET,0. ,0.0, I ,XP(2) ) 
(-TS,ET,T5.n. ,9.0,1,»P(3) ) 
(-TS,ET ,TS,n. ,ES2. 1 , XP (A) > 
(-TS,ET,TS,0.,«.0,2,xP(5)l 
(-TS,ET,TS,0.,£S2,2,xP(6I) 
(‘TS,0. ,0. ,8. . 5. 0, I , XP(7)  ) 
(♦TS,0. ,0. ,8. ,ES2, 1 ,XP (8) ) 


PI  = XP(1)  ♦ (XP(3)  ♦ xP(5))/2.  ♦ XPin 
PI  = PI  - XP(2)  - (XP(A)  ♦ XP(A)l/2.  - XP(8) 
35  PRINT  37,  XMU ,EB, TB , ES , T S 
37  FORMAT  (1X.5E20.7//) 


PRINT 

AO 

, XP ( 1 ) , XP ( 3)  , 

XP(5) , 

XP(7) 

FORMAT  (1X.AE20.7/) 

PRINT 

40 

, XP(2),  XP(A), 

XP(61 , 

XP(8) 

XP(I) 

= 

XP(1)  - XP(2) 

XP(3) 

= 

( XP(3)  - XP(A) 

) / 2. 

XP(5) 

= 

( XP(5)  - XP(6) 

) / 2. 

XP(7) 

= 

XP(7)  - XP(8) 

PRINT 

50 

, XP(1) , XP(3) . 

XP(5) . 

XP  (7) 

50  FORMAT  (IX,5E20.7////) 
55  CONT  IN  IE 
STOP 
END 


60  SUBROUTINE  EYEL  ( A , B , C , 0 • U , L , X I ) 
DIMENSION  xx(3l 
COMMON  XK 
CALL  BASIC  (C,0,U) 

XI  = XK(L1 

call  BASIC  (A,B,U) 

XI  = XI  - XK(L) 

RETURN 

END 


C 

350  SUBROUTINE  BASIC  (G,H,U) 

DIMENSION  XK(3) 

COMMON  XK,  TS 

X2  = SORT  ( I . . G»G) 

X3  = (G»H)/X2 
Y = X2»U  ♦ X3 

YN  = (.a86??693)  • EXP(-J*U)  • FRF(G“U  ♦ M) 
XA  = EXP  (- (H»H)  / (X2»X2)  ) '' (2,*X2) 

X5  = (1.772A539)®EDF(y) 

XMO  = XA  • X5 
IF  (Y.LE.12.)  go  to  355 


XMl  = -(X3  • XM0I/X2 

XM2  = (X5«  ( I .♦2.»X3»X3' •»<•)  / (2.»X2»X2) 

GO  TO  360 

355  XMl  = - (XA»FXP  (-y»y)  . x3*XM0)/X2 
XM2  = 2.®(?.®X3  - Y)  • f*^(*Y”Y) 

XM2  = ((Xm?  ♦ X5®(1.  • 2,*X3«X3I )•XA)/(^.•X2•X2I 
360  XK(1)  = (I. 1283792)  • (G*«XMO  ♦ XM2)  - YN  • (1.  ♦ U*U) ) 
XK(2)  = -(XM2  • .56A  1 8RX,)  / rs 

XK(3)  = ( ( .56A1896) / ( TS® >S) ) • (G*XMO  - YN  - 2.»(H»XM1  , 

return 

END 

//G.SYSIN  On  • 

♦ 2, 666666 7 •OA. 0 00 ♦OH. 000 
/• 


1 

t 


f.»XM2)  ) 
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C TWO  COLLISIONS.  XMu  INPUT  FOR  4/3,  8/J 

C 

DIMENSION  *K(3),  XP(lft),  P20(3).  PI(3),  QI(3),  RI(3).  XKI(41 
COMMON  XK,  TS 
C INPUT 

DO  6b  J = 1 .26 

4 READ  200.  XMU.EB.TB.EPR 

£S  = SORT  ( (I.  ♦ XMU  » • FB  ) 

ES2  = ES/2. 

TS  = 1.  / S0RT(XMU»TB) 

ET  = ES  • TS 

EB2  = SORT  (FBR) 

UA  r ( (3. -XMU)  / (1.  ♦ XMU)  ♦ EB2  ) • ES  • .5 

V2R  = ( (3.  • XMU  - 1.)  • TS)  / (3.  - XMU) 

V2I  = -(  (1.  ♦ XMU)  • Eb  • TS  « E82)  / (3.  - XMU) 

P2B(3)=  0. 

PRINT  20*.,  XMU,  EB.  TB.  ES,  TS.  EBR 
PRINT  200 
C 

5 Z1  = ( XMU  - 1 . ) • TS  • .S 
IF  (EBR.LE.l.)  GO  TO  10 


uc 

= (2 

. • ES  • EB2) 

/ (1. 

♦ XMU) 

CALL 

EYEL 

(-TS, 

ET,  V?R, 

V2I , 

UC. 

1 , 

XP(  1 ) 

CALL 

ETEL 

(-TS. 

ET,  V?P, 

V2I  , 

UA, 

1 . 

XP(2) 

CALL 

ETEL 

(-TS. 

ET,  V?R, 

V2I  , 

UC. 

2. 

XP  ( 3) 

CALL 

ETFL 

(-TS, 

ET,  V?R. 

V2I, 

UA . 

2. 

XP  (,,  ) 

CALL 

EYEL 

(-TS. 

ET.  21  . 

0.0. 

B.  . 

1 . 

xp(b) 

CALL 

EYEL 

(-TS, 

ET,  2!  , 

0.0, 

UC. 

1 . 

XP  (6) 

CALL 

EYEL 

t-TS. 

FT.  21  . 

0.0, 

B.  . 

2, 

XP(  7) 

CALL 

EYEL 

(-TS, 

ET,  21  , 

0.0, 

UC, 

2. 

XP(H) 

PRINT  206,  XP(1),  XP(3),  XP(5),  XP(7) 

PRINT  206,  XP(2),  XP(4),  XP(6),  XP(B) 

DO  6 K 1 = 1.7,2 
6 XP (K1 ) = XP(K 1 ) - XP  (K1  . 1 ) 

PRINI  206,  XP(1),  XP(3).  XP(5).  XP(7) 

P2A  = ( XP(I)  ♦ XP(3)  ♦ XP(5)  ♦ XP(T)  ) • .5 
PRINT  210,  P2A 
GO  TO  15 


22 

= (2 

.•ES) 

/ (1.*xmJ) 

CALL 

EYEL 

(-TS, 

ET.  21  , 

0. 

f B • t 

1 , 

XP(  1 ) 

CALL 

EYEL 

(-TS, 

ET,  21  . 

0. 

. 7?. 

I , 

XP(2) 

CALL 

EYEL 

(-TS. 

ET,  21  , 

0. 

f B • • 

2. 

XP(3) 

CALL 

EYEL 

(-TS, 

ET.  21  , 

0. 

. 72. 

2, 

XP  (4  ) 

PRINT 

212, 

XP(1) 

, XP(3) 

print 

212, 

XP  (2) 

, XP(4) 

XP(1)  = XP(1)  - XP(2) 

XP(3)  = XP(3)  - XP(4) 

PRINT  212,  XP ( 1 ) , XP ( 3) 

P2A  = ( XP(1)  . XP(3)  ) • .5 
PRINT  214,  P2A 

15  IF  (XmU.GE.3.)  go  to  f>5 
IF  (XMU.LE.  1 .667)  GO  TO  <?5 


20  PI ( 1 ) = 

1 ./4. 

PK2)  = 

0. 

PK3)  = 

0. 

OKI)  = 

-3./1 1 

OK2)  = 

-3./1 1 

01(3)  = 

0. 

RI  ( 1 ) = 

21 ./44 

RK2)  • 

0./1 1 . 

RIO)  = 

0. 

XKI  ( 1 ) = 

S.  /6. 

XKl (2) = 

1.0 

XKI (J)= 

0./3. 

XKl (4)1 

0. 

IL 

2 

GO  TO  2« 


47 


c 


c 


c 


c 


25  PI  (1)  = I./**. 

PI(2)  * -5./?e. 

PI(3)  = 0. 

OKI)  = -3./7. 

QI (2)  = 

01 (3>  = - 3./7. 

RI  ( 1 > = 9./P8. 

RK2)  = 0. 

RIO)  = 8./7. 

XKI(1)=  1./6. 

XKI(2)=  3./8. 

XKI<3)=  l.n 
XKI(8)=  4./3. 

IL  =3 

28  PRINT  216t  PI 
PRINT  216.  01 
PRINT  216,  RI 
PRINT  206,  XXI 
PRINT  208 

DO  60  I = 1 ,1L 

EBRI  = ( (2./(l.*XKI (I) ) ) - ( (3,-XMJ)/(l.»XNU) ) ) ••  2. 

EBRIl  = ( (2./(l.*XXI (1,1) ) ) - { (3.-XMU)/(l.«XMU) ) ) *•  2. 

TSK  = TS  • XKl  ( I ) 

TSKl  = TS  * XXI  (I»l) 

ESK  = ES  / (l.»XXl  (1)  ) 

ESKI  = ES  / (1.  ♦ XX1(1,1)  ) 

UC  = ( (1.,XM0)«ES»ER2)  / ( (3.«XMU  - 1.)  - ( 3.-XMU) «XxI  ( I«  1 ) ) 

UE  = ( ( 1 .‘XMU) ‘ES'ESa) / ( <3.«XMU-1 . ) -(3,-XMU) »XXI  ( I ) ) 

PRINT  216,  EBRIl,  EBR,  E8R1 
PRINT  212,  UE.UC 

30  IE  (EBR. GE. EBRI 1 ) GO  TO  35 


CALL 

EYEL 

(TSX. 

0.  . 

TSXl , 

0.. 

8.0 

3, 

XP(1)  ) 

CALI 

EYEL 

(TSX, 

0., 

TSXl , 

0., 

ESX 

3. 

XP(2)  ) 

call 

EYEL 

(TSX, 

0.  . 

TSXl  . 

0., 

8.0 

2, 

XP(3)  ) 

CALL 

EYEL 

(TSX, 

n.. 

TSXl, 

0., 

ESK 

2, 

XP(8)  I 

CALL 

EYEL 

(TSX, 

0., 

TSXl  , 

0 • « 

8.0 

1, 

XP(5)  ) 

CALL 

EYEL 

(TSX, 

0. , 

TSXl  , 

0 • « 

FSX 

1 , 

XP(6)  ) 

CALL 

EYEL 

(-TS, 

ET, 

TSXl  . 

0., 

ESX 

3, 

XP(7)  I 

CALL 

EYEL 

(-TS, 

ET. 

TSXl, 

0.  . 

FSXl 

3. 

XP(8)  ) 

CALL 

EYEL 

(-TS, 

ET. 

TSXl  , 

0. , 

ESK 

2, 

XP(9)  ) 

CALL 

EYEL 

(-TS, 

ET, 

TSXl , 

0., 

ESXl 

2, 

XP(IO)  ) 

CALL 

EYEL 

(-TS, 

ET, 

TSXl, 

0 • « 

rSx 

1 , 

XP(ll)  ) 

CALL 

EYEL 

(-TS, 

ET, 

TSXl  , 

0., 

ESX  1 

1, 

XP(12)  ) 

PRINT  218 
GO  TO  38 


35  IE  (EBR. LE. EBRI ) GO  TO  85 


CALL 

EYEL 

(TSX, 

0., 

TSXl  , 

0.0, 

8., 

3, 

XP(  1 ) 

CALL 

EYEL 

(TSK, 

0. . 

TSXl  , 

0.0. 

UC. 

3, 

XP(2) 

CALL 

EYEL 

(TSX, 

0 • « 

TSXl  , 

0.0. 

8., 

2, 

>P(3l 

CALL 

EYEL 

(TSX, 

0 • • 

TSXl. 

0.0, 

UC. 

2, 

XP(8) 

CALL 

EYEL 

(TSX, 

0. . 

TSXl , 

0.0, 

8.. 

!• 

XP(5) 

CALL 

EYEL 

(TSX. 

0.  , 

TSXl  , 

0.0, 

UC, 

1 , 

XP  (6) 

CALL 

EYEL 

(TSX. 

0 • • 

V2R  , 

V2I. 

UC, 

3, 

XP  ( 7) 

• CALL 

EYEL 

(TSX, 

n. , 

V2R  • 

Y2I  . 

UE, 

3, 

XP(8) 

CALL 

EYEL 

(TSX, 

0. . 

V2R  , 

V2I. 

UC. 

2, 

XP  (9) 

CALL 

EYEL 

(TSX, 

0.  . 

V2P  , 

V2I. 

UE. 

2. 

XP ( 1 0) 

CALL 

EYEL 

(TSX, 

n.. 

V2R  • 

V2I  , 

UC. 

1, 

XP( 1 1 I 

CALL 

EYEL 

(TSX, 

0.  . 

V2»  • 

Y2I, 

UE  , 

1* 

XP( 12) 

PRINT  2?0 

38  PRINT  208.  XP(1),  XP(3).  XP(S),  XP(7).  XP(9).  XP(1H 

PRINT  208,  XP(2).  XP(8),  XP(6),  XP(8),  XP(IO),  XP(12) 

DO  80  XI  * 1.11.2 

80  XP(Xl)  I XP(XI)  - XP(Xl.l) 

PRINT  208,  XP(1),  XP(3).  XP(5),  XP(7),  XP(9),  XP  ( 1 II 

DO  82  XI  = 1,7,6 

XP(XI I = XP(X 1 I • PI ( 1 1 

XP(Xl ,2)  r XP (XI »2l  • 01 ( I ) 

82  XP(Xl«8)  = XP(Xl,8l  • RI(II 
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PaB(l)  = XP(l)  ♦ XP(3)  ♦ XP(5)  » XP(7)  ♦ XP(9)  • XP(U) 
PRINT  ?26,  P2B(I) 

GO  TO  60 
C 

<*5  EBPR  = ( (3.«XMU-1 . )- (3.-XBU)»XKI  ( ) 

E8RR  = I EBRP  / I (1.*XMJ)*(1.*XKI (I) ) ) I ••  2. 

PRINT  222,  EBRR 


IE  (EBR.GE.E0RR)  GO 

TO  50 

CALL 

ETEL 

(TSK, 

0., 

TSKl , 

0.0, 

0.0, 

3, 

XP(1) 

CALL 

EYFL 

(TSK. 

0 , , 

TSKl  , 

0.0, 

ESK, 

3, 

XP(2) 

CALL 

EYEL 

(TSK, 

0 . , 

TSKl  , 

0.0, 

0.0, 

2, 

XP13) 

CALL 

EYEL 

(TSK, 

0.  , 

TSKl  , 

0.0, 

ESK , 

2, 

XP(<,) 

CALL 

EYEL 

(TSK, 

0.  , 

TSKl  , 

0.0, 

8.0, 

1 , 

XP(5) 

CALL 

EYEL 

(TSK, 

0.  , 

TSKl  , 

0.0, 

E SK  , 

1 , 

XP(6) 

) J 

CALL 

EYEL 

(-TS, 

ET, 

TSKl  , 

0.0, 

FSK  , 

3, 

XP(  7) 

CALL 

EYEL 

(-T5, 

ET. 

TSKl  , 

0.0, 

UC  , 

3, 

XP(8) 

CALL 

EYEL 

(-TS, 

ET, 

TSKl  , 

0.0, 

ESK, 

2. 

XP(9) 

CALL 

EYEL 

(-TS, 

ET, 

TSKl  , 

0.0, 

UC  , 

2, 

XP( 10) 

) 

CALL 

EYEL 

(-TS, 

ET, 

TSKl  , 

0.0, 

ESK. 

1 , 

XP(ll) 

) 

CALL 

EYEL 

(-TS, 

ET, 

TSKl  , 

0.0, 

UC  , 

1 , 

XP(12) 

) 

CALL 

EYEL 

(-TS, 

ET, 

V?R  , 

V2I. 

UC  , 

3, 

XP( 13) 

) 

CALL 

EYEL 

(-TS, 

ET, 

V2R  , 

Y2I, 

UA  , 

3, 

XP(IS) 

) 

CALL 

EYEL 

(-TS, 

ET, 

V2R  , 

Y2I, 

UC  , 

2, 

XP (15) 

) 

CALL 

EYEL 

(-TS, 

ET, 

V?R  , 

V2I, 

UA  , 

2, 

XP ( 16) 

) 

CALL 

EYEL 

(-TS, 

ET, 

V2R  , 

Y2I, 

UC  , 

1 , 

XP( 17) 

) 

CALL 

EYEL 

(-TS, 

ET, 

V2R  , 

V2I, 

UA  , 

1, 

XP( 18) 

) 

GO  TO  55 

50  CALL 

EYEL 

(TSK, 

0. 

TSKl , 

0.0, 

0.0, 

3, 

XP(  1 ) 

) 

CALL 

EYEL 

(TSK, 

0. 

TSKl  , 

0.0, 

UC  , 

3, 

XP(2) 

) 

CALL 

EYEL 

(TSK, 

0. 

TSKl  , 

0.0, 

8.0, 

2, 

XP(3) 

) 

CALL 

EYEl 

(TSK, 

0. 

TSKl , 

0.0, 

UC  , 

2, 

XP(4) 

) 

CALL 

EYEL 

(TSK, 

0. 

TSK  1 , 

0.0, 

8.0, 

1, 

XP(5) 

) 

CALL 

EYEL 

(TSK, 

0. 

TSKl, 

0.0. 

UC  , 

1, 

XP(6) 

) 

CALL 

EYEL 

(TSK, 

0. 

V2R  • 

V2I. 

UC  , 

3, 

XP(7) 

) 

CALL 

EYEL 

(TSK, 

0. 

V2R  , 

V2I, 

ESK, 

3, 

XP(8) 

) 

CALL 

EYEL 

(TSK, 

0. 

V2R  , 

V2I, 

UC  , 

2, 

XP(9) 

1 

CALL 

EYEL 

(TSK, 

0. 

V2R  , 

V2I, 

ESK, 

2. 

XP(IO)  ) 

CALL 

EYEL 

(TSK, 

0. 

V2P  , 

V2I. 

UC  , 

1, 

XP(ll)  ) 

CALL 

EYEL 

(TSK, 

0. 

V2R  , 

V2I , 

ESK, 

1, 

XP(12)  ) 

CALL 

EYEL 

(-TS, 

ET, 

V2R  , 

V2I, 

ESK, 

3, 

XP( 13) 

) 

CALL 

EYEL 

(-TS, 

ET, 

V?R  , 

V2I, 

UA  , 

3, 

XP ( 19) 

CALL 

EYEL 

(-TS, 

ET, 

V?R  , 

V2I. 

FSK, 

2, 

XP ( 15) 

) 

CALL 

EYEL 

(-TS, 

ET, 

V?P  , 

V21, 

UA  , 

2, 

XP(16) 

CALL 

EYEL 

(-TS, 

ET, 

V?R  , 

V2I, 

FSK, 

1 , 

XP( 1 7) 

CALL 

EYEL 

(-TS, 

ET, 

V2R  , 

V2I  , 

UA  , 

1 , 

XP( 18) 

) 

S5  PRINT  209,  XP(1),  XP(3),  XP(5).  XP(7),  XP(9),  XPIll) 

PRINT  209,  XP(2),  XP(i.»,  XP(6),  XP(8),  XP(IO),  XP(I2I 
PRINT  216,  XP(13),  XPtlS),  XP(17> 

PRINT  216,  XP(19),  XP(16»»  XP(I8) 

PRINT  200 
DO  57  K1  = 1 ,17,2 
57  XP(Kl)  = XP(Kl)  - XP(Kl,l) 

PRINT  209,  XP(1>.  XP(3»,  XP(5».  XP(71,  XP(9>,  XP ( 11 ) J 

PRINT  216,  XP(I3),  XP(15),  XP(l7) 

DO  59  K 1 = 1,13,6 
XPIKl ) = XP(Kl)  • PI  ( n 
XP(K1*2)  = XP(K],2)  » OKI) 

59  XP(Kl,9)  = XP(Kl*9)  • P1(II 

P2BII)  = XP(1)  ♦ XP(3)  ♦ XP(5)  ♦ XP(7»  ♦ XP(9) 

P2B(I)  = XP(ll)  ♦ XP(13)  ♦ XP(15)  ♦ XP(17)  ♦ P2R ( I) 

PRINT  226,  P2B(1) 

C 

60  CONTINUE 

P2BB  = P2a  ♦ P20ni  ♦ P2H(2)  ♦ P2H(3) 
print  2?9,  P2PR 
65  CONTINUE 
STOP 
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200 

204 

206 

zoe 

210 

212 

214 

216 

218 

220 

222 

224 

226 


C 

300 


C 

350 


format 

format 

format 

format 

format 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

END 


(Fln.7.F6.3.E8.1,E11.4l 
(lx,  6E20.7) 

(lx,  4E20.7) 

(lx,//) 

(IX,  »P2A  EBR  GREATER  THAN  ONE  =*t  £20. 7//) 

(IX,  2E20.7) 

(lx,  »P2A  EBR  LESS  THAN  ONE  =••  E20.7//) 

(lx,  3E20.7) 

(IX,  *P29  EBR  LESS  THAN  E0RI1" 

(lx,  *P2B  EBR  greater  Than  EBRIm 

(IX,  *P?B  EBR  bracketed  EBRR  ='.  E 20.7  ) 

(lx.  *TW0  COLLISION  PROBABILITT  = •,  E 20.7/////) 
(IX,  •P2BI  = ',  £20.7//) 


SUBROUTINE  EYEL  (A,B,C«0*U,L,XI ) 
DIMENSION  *K(3) 

common  XK 

CALL  BASIC  (C,0,U) 

XI  = XK(L) 

CALL  BASIC  (A.B.U) 

XI  = XI  - XK(L) 

RETURN 

END 


SUBROUTINE  BASIC  (G.H.U) 
DIMENSION  XK(3) 

COMMON  XK,  TS 

X2  = SORT  (1.  ♦ G*G) 

X3  » (G»H)/X2 


Y = X2*U  ♦ X3 

YN  = (.886?2693>  • EXP(-U»U)  • ERF(6»U  ♦ H) 

X4  = EXP(-(H»H)/(X2«X2) )212.»X2) 

X5  = (1.772453R)»ERE(Y) 

XMO  = X4  * X5 
IF  (Y.LE.12.1  go  to  355 
XMl  = -(X3  • XM0)/X2 

XM2  = (X5»(1.*2.«X3*X3>*X4)  / (2.»X2*X2) 

GO  TO  360 

355  XMl  = - (X4»EXP(-Y«Y)  • x3»XM0)/X2 
XM2  = 2.»(2.»X3  - Y)  * EXP(-Y*Y) 

XM2  = ((XM?  ♦ X5»(l.  ♦ 2.»X3*X3))«X4) /(2,»X2*X2) 

360  XK(1)  = (1.1283792)  • (G'fXMO  ♦ XM2)  - YN  • (1.  ♦ U»U)) 
XK(2)  = -(XM2  • .5641B98)/TS 

XK(3)  = ( (.5641896)/(TS»TS) ) • (G»XMJ  - YN  - 2.»(H»XM1  . 

RETURN 

END 

//G.SYSIN  OD  • 

♦ 1. 3333333* 0.000 *5. OE-01*O.0n00E» 00 

/• 
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SAMPLE  OUTPUT.  TWO  COLLISION  PROGRAM 
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TABLE  2 

ONE  COLLISION  PROBABILITY,  y = 4/3 


T 

f’l 

e-|  = 0.5 

1 .0 

2.0 

3,0 

4.0 

5.0 

0.25 

.693 

.405 

.110 

.0261 

.00575 

.00121 

0.50 

.711 

.463 

.168 

.0548 

.0168 

.00495 

1.00 

.741 

.541 

.267 

.123 

.0548 

.0237 

2.00 

.783 

.630 

.397 

.243 

.146 

.9863 

4.00 

.829 

.716 

.533 

.393 

.287 

.209 

6.00 

.854 

.761 

.605 

.480 

.380 

.299 
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TABLE  3 

TWO  COLLISION  PROBABILITY,  y = 4/3 


T 

Pll(ei, 

0) 

Si  = 0 

1 

1 .5 

2 

3 

4 

0.1 

.872 

.307 

.145 

.0644 

.0115 

.00188 

0.25 

.791 

.286 

.147 

.0726 

.0164 

.00346 

0.5 

.715 

.265 

.150 

.0828 

.0244 

.00697 

1.0 

.623 

.239 

.148 

.0920 

.0357 

.0140 

2.0 

.517 

.207 

.138 

.0941 

.0453 

.0226 

4.0 

.409 

.170 

.120 

.0864 

.0476 

.0273 

6.0 

.350 

.148 

.106 

.0784 

.0453 

.0274 

T 

Pii(ei  , 

1) 

PllOO'^ 

1 ^ 

1 u ! 

Cl  - 1 

2 

3 

4 

£2  “ 1 

2 

0.1 

.288 

.0640 

.0115 

.00188 

.544 

.225 

0.25 

.241 

.0669 

.0158 

.00342 

.503 

.222 

0.5 

1 

.201 

.0660 

.0205 

.00620 

.460 

.217 

1.0 

.162 

.0606 

.0235 

.00934 

.400 

.202 

2.0 

.125 

.0512 

.0229 

.01067 

1 

.328 

.175 

4.0 

.0937 

.0404 

.0195 

.00991 

.255 

.141 

6.0 

.0781 

.0343 

.0170 

.00891 

.215 

.120 

54 


TABLE  4 

ONE  COLLISION  PROBABILITY,  p = 8/3 


T 

= 4.0 

8.0 

16. 

32. 

.5 

.00780 

.000026 

.000000 

.000000 

1.0 

.0431 

.00121 

.000000 

. 000000 

2.0 

.137 

.0167 

.000217 

.000000 

4.0 

.284 

.0833 

.00737 

.000058 

6.0 

.379 

.154 

.0275 

.000938 

8.0 

.445 

.216 

.0553 

.004017 
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TABLE  5 

TWO  COLLISION  PROBABILITY,  y = 8/3 

T I 


0 

II 

1 

3 

6 

8 

.910 

.00565 

.00032 

.00002 

.833 

.0230 

.00353 

.00053 

.732 

.0542 

.0156 

.00450 

.613 

.129 

.0329 

.0135 

.539 

.133 

.0403 

.0186 

.488 

.130 

.0432 

.0212 

Pll(ei. 

£2) 

F-,  = 10’^ 

4 

4 

4 

8 

8 

£2=4 

2 

4 

8 

4 

8 

.0889 

.00084 

.0875 

.00529 

.00344 

.000455 

.000016 

.000006 

.0834 

.01872 

.00833 

.000629 

.000312 

.000047 

.0756 

.0399 

.0136 

.000734 

.00180 

.000137 

.0652 

.0581 

.0172 

.000757 

.00457 

.000237 

,0537 

.0641 

.0181 

.000708 

.00694 

.000297 

.0388 

.0553 

.0156 

.000572 

.00747 

.000294 

.0290 

.0438 

.0124 

, 000450 

.00641 

.000250 
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